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Abstract: Temperature measurement with nano-Kelvin resolution is 
demonstrated at room temperature, based on the thermal dependence of an 
optical crystal anisotropy in a high quality whispering gallery resonator As 
the resonator's TE and TM modes frequencies have different temperature 
coefficients, their differential shift provides a sensitive measurement of 
the temperature variation, which is used for active stabilization of the 
temperature. 
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1. Introduction 

Whispering Gallery Mode (WGM) resonators hold a unique place among optical resonators 
for several reasons. Relying on continuous total internal reflections for light confinement, such 
resonators feature extremely high quality factors Q{Q> lO'^ has been demonstrated (T||) in a 
very broad optical range. They are also very compact, with a typical size ranging from a few 
tens of microns to several millimeters, and are entirely solid-state. It is therefore not surprising 
that among many other applications, WGM resonators have been suggested for use as reference 
cavities for laser stabilization E |3] |4] 2j l6i . Such an approach to laser stabilization is predicted 
IFTI to allow for the fractional optical frequency stability at the level of one part per 10^'^ at 1 
s. In addition to optical clocks, laser stabilization with a WGM resonator holds a great promise 
for high-resolution spectroscopy, optical metrology and other fields. 

The main practical difficulty associated with using WGM resonators for laser stabilization 
arises from the fundamental fact that the light propagates inside an optical material. This ap- 
pears to be a disadvantage compared with e.g. low-expansion Fabri-Perot resonators that are 
essentially vacuum-filled. In WGM resonators the material's optical properties set the limit not 
only for the resonator's g-f^ctor, but also for its stability. Indeed, it is easy to see that even a 
very small variation A(nL) of the resonator's optical length nL such that A{nL)/{nL) = 1/2 
causes the resonance frequency to shift by one linewidth. 

The most important factors contributing to variation of A(«L) are the thermal refractivity and 
thermal expansion |7J. One approach proposed iHH] for compensating these variations is to 
use the temperature dependence of the resonator's own anisotropy. In this work we report the 
first practical implementation of this approach and demonstrate active temperature stabilization 
of a WGM resonator at 76.04 °C. We achieve temperature stabilization at the level of 200 nK 
when integrated for 1 second, and below 10 nK when integrated for 10,000 seconds. We see 
this result as significant not only in the context of laser stabilization but also as a demonstration 
of a novel temperature sensor with a broad spectrum of potential applications. 

2. Theoretical estimates 

Temperature dependence of a WGM frequency / can be found by differentiating the WGM 
dispersion equation (found in e.g. |9!l) with respect to temperature T: 

1 d f I dn I dR _ 
fdf^ndf^Rdf^^' 

In Eq. ([T]i is the resonator radius and n is the refractive index for the given polarization and 
wavelength X. This equation is an approximation, where the small terms of the order X /R and 
higher have been neglected. 

The second and third terms of Eq. ([T) are equal to the thermorefractivity and thermal ex- 
pansion coefficients, aj,"'''^ and aj"^ respectively. Here the superscripts indicate ordinarily and 
extraordinarily polarized light. Since our resonator is shaped as a disk made from a z-cut mag- 
nesium fluoride (MgF2) wafer, its TE and TM modes support the ordinary (horizontally) and 
extraordinary (vertically) polarized fight, respectively. Therefore the thermorefractivity of the 
former is approximately equal to oci"', and of the latter, to a^^\ The thermal expansion coeffi- 
cient, on the other hand, is a^"' for both mode families. These coefficients for MgF2 are given 
by Q: 

ajf'^ ^ (7.14-0.040627) X 10"^ 



a,r = (3.02 - 0.04071 r) X 10"^/:"^ 

a/"' = 9xlO-''K-\ (2) 

where T is in °C. From we see that while the thermal variation of each mode's frequency is 
dominated by the thermal expansion, 

■^fo.e « - a/"' « - 1 .73 GHz/K, (3) 

the variation of the frequency difference Af ~ f,, — fe depends only on the differential ther- 
morefractivity; 

J; A/ = -j (ai") - ai^)) « -79. 1 MHz/K. (4) 
We use this temperature dependence to measure and stabilize the resonator temperature. 

3. Experimental demonstration 

Our experimental setup is shown in Fig. [1] The key element of this setup is a MgF2 WGM 
resonator which is shaped as a disk of approximately 8 mm in diameter and 2 mm thick. Its 
modes are excited by a CW DFB laser at A = 1560 nm. Using a fiber optic polarization adjuster 
we set input light polarization such as to excite both TE and TM mode families. The optical 
power at the output of the polarization adjuster ranges from 0.5 to 1.5 mW in each polarization. 

The light is coupled into and out of the resonator using the standard /rMifrafec/ total internal 
reflection technique with an angle-polished fiber |10| . The coupling rate is controlled using a 
micro screw and a two-springs differential that can adjust the distance between the resonator 
rim and the fiber coupler to a fraction of the optical wavelength. 

The resonator and fiber coupler are mounted inside a 1 inch brass cube, shown in the photo 
in Fig.[T] This cube also has two built-in resistive heaters and a thermistor to provide the initial 
temperature stabilization by means of a commercial temperature controller During the active 
phase of our experiment this controller is disengaged and the heaters are controlled by a feed- 
back loop described below. The cube is placed in a cardboard box on an optical table, with no 
other means of external temperature stabilization. 




Fig. 1. The experimental setup. Temperature of the brass cube on the photo is set to 76 °C 
and stabilized to a few hundred nanodegrees using the dual-mode stabilization technique. 




Fig. 2. (a) A pair of TE and TM WGMs withi a small frequency detuning A/ = /o — /c- 
(b) Temperature dependence of Af{T) is linear with a -89.8 MHz/K slope. The inset shows 
a small distortion around the set point Af{Tsei) = which is due to the WGMs cross- 
coupling. 



The out-coupled light is coUimated by a lens and split by a polarizing beam splitter. Two 
photodetectors at two ports of the beam splitter can therefore observe the TE and TM modes. 
The photodetectors' signals are acquired by a computer-controlled DAQ which also controls 
the laser frequency and performs the temperature stabilization. 

By varying the laser wavelength or the resonator temperature it is possible to find multiple 
pairs of TE and TM modes that have sufficiently close frequencies /„ and f^, as shown in 
Fig. 12 a). In our resonator a typical intrinsic (weakly coupled) value of Q for TE WGM is 
2 X 10^ and for TM WGM is 1 x lO'*. 

To observe the WGMs, the DAQ sweeps the laser frequency 100 times per second, each time 
acquiring 2500 data points per channel. This data is analyzed in real time during each sweep, 
and the WGM frequencies /„ and fe are found by peak fitting. Their sum is used to continuously 
re-center the laser frequency sweep with respect to the modes. Their difference A/ = — 
fe is used as the temperature error signal. This signal has been calibrated by measuring Af 
simultaneously with the resistance of the thermistor, which required a seven-digit Ohm meter 
The calibration curve shown in Fig. |2lb) yields the temperature tuning rate of d{Af)/dT = 
—89.8 ± 0.2 MHz/K, which is reasonably close to the theoretical value found in 

To demonstrate the use of our temperature measurement technique for temperature stabiliza- 
tion, we have implemented two Labview-based feedback loops with proportional and integral 
gains. The slower loop uses the heaters mounted inside the brass cube to provide the long-term 
stabilization. The faster loop uses an amplitude electro-optical modulator to vary the optical 
power, which has a high-speed heating effect ifTTl lTS 13 1. The proportional and integral gains 
of both feedback loops, as well as the cross-over time of the loops, are optimized for the best 
short term and long term performance. 

When the temperature control loops are closed, the TE and TM modes shown in Fig. |2a) 
are brought on resonance Af{Tset) = as shown in Fig.[3ta), and the resonator temperature is 
stabilized at that point. Special care needs to be taken to choose the WGMs that have weak 
cross-coupling, which may arise from either linear or nonlinear interaction between the WGMs 
and can severely distort the WGM spectrum (see Fig. lUb)). Linear cross-coupling may be 
caused by inhomogeneities, such as e.g. surface scratches or dust. Such defects are not usually 
expected in high-Q resonators. Nonlinear coupling may be due to thermal, acoustical or Kerr 
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Fig. 3. (a) The WGMs shown in Fig. |2]are brought on-resonance by temperature control. 
These modes have weak nonlinear coupling, (b) Another pair of modes couples strongly 
and is not suitable for our application. 



nonlinearities, and is directly related to the WGMs spatial overlap. 

The modes shown in Fig. |3ja) interact only weakly. To verify that this interaction does not 
significantly alter the linear dependence Af{T) near the set point, we took a number of data 
points within the WGM linewidth. These data points are shown in the inset in Fig.|2jb). With 
such a fine step we should be able to capture all reasonably expected distortions. We see that 
while the cross-coupling effect is clearly present, it is not likely to strongly affect our measure- 
ments. 

To characterize the performance of our temperature stabilization setup, we record a time 
series of the Af fluctuations and find both their spectral density and Allan variance. These are 
converted to the temperature units using the experimental value for d{Af) /dT . A typical result 
is shown in Fig. |4] The optical powers injected in the WGM resonator in this case were 0.8 
mW for TE and 1.1 mW for TM polarizations. The peculiar non-monotonic behavior of the 
Allan variance below 5 seconds that can be seen in this Figure is a consequence of a very high 
integral gain set in the slow control loop, which made it nearly unstable. The high gain was 
necessary to provide a good long-term stability against the ambient temperature variations to 
which our system has been considerably exposed. The residual signature of these variations is 
visible around 1000 seconds and is likely to be related to the building air conditioner, which 
has an appropriate period of operation. Apart from this signature, the Allan variance in Fig. |4] 
scales as 0.5 /xK / s'/^ for the averaging time longer than 5 seconds. The inverse square-root 
behavior of the Allan variance has been theoretically predicted in Q . 

For comparison, we measured the temperature Allan variance of the resonator stabilized by 
the commercial digital temperature controller, and found it roughly flat at a level of 0.2 to 0.7 
mK. Therefore, the temperature stabilization technique we report outperforms this controller by 
over three orders of magnitude with only 1 second integration time. It is also interesting to com- 
pare our results on temperature stabilization with sensitivity of the state-of-art non-cryogenic 
temperature sensors (which are sensing but not controlling temperature) developed for LISA 
lfT4l . This sensitivity is given as the noise-equivalent temperature spectral density, which is 
also shown in the inset of Fig.|4] From this we see that our temperature stabilization technique 
considerably outperforms the LISA sensors. 

The temperature sensitivity of our approach is apparently limited by the accuracy of deter- 
mining the WGM resonance frequencies. A least-square Lorentzian fitting of the peaks such as 




Fig. 4. The Allan variance of the temperature fluctuations in the stabilized WGM resonator 
in nano-Kelvin and in terms of the expected fractional frequency stability. The straight line 
has a slope of 0.5 flK / s'/-. The error bars are based on the measurements statistics. On the 
inset: the spectral density of the temperature fluctuations measured in our system (below) 
and the noise-equivalent temperature spectral density in 1 141 (above). 



shown in Fig. [3ta) yields the standard error of their frequency difference which corresponds to 
OA/(10ms) 4.8 /iK. With 100 such measurement per second we may expect C7a/(1 s) ~ 480 
nK. The observed one-second Allan variance of less than a half of this value indicates that we 
are at or near the accuracy limit of our measurement. It should be pointed out that the WGM 
peaks in Fig. |3la) are noticeably asymetric due to the dynamic heating by the laser sweep 
ifTTl [T2I [T3I . Therefore reducing the sweep time with a faster DAQ would not only increase 
the data rate but also improve each peak fitting accuracy by suppressing this asymetry, and will 
improve the temperature sensitivity of our method. 

Since the main motivation of the demonstrated temperature stabilization technique is to pro- 
vide a stable laser locking reference, it is interesting to evaluate the projected fractional fre- 
quency stability of such a reference. This can be done by multiplying the measured temperature 
stability by the coefficient (|3j. This result is also shown in Fig.|4] According to this, we should 
be able to stabilize our laser to one part per 7 x 10^'^ per 1 s integration time, and down to 
one part per 6 x 10^^^ per 10,000 s integration time even with the present proof-of -principle 
system. 

It should be pointed out that in addition to its temperature dependence, the dual-mode fre- 
quency detuning A/ may have other dependencies, e.g. on the atmospheric pressure, TE/TM 
power ratio, etc. It also could be time-dependent due to crystal aging. Successfully using A/ 
for temperature stabilization however shows that its temperature dependence is dominant, and 
serves as a validation of our technique of ultra-sensitive temperature measurements. A further 
verification of the achieved temperature stability will be available from measuring the rela- 
tive frequency stability of two independently stabilized lasers, which will be the subject of our 
future work. 



4. Conclusions 



We have demonstrated a WGM based temperature sensor with nano-Kelvin sensitivity oper- 
ating at room temperature. We have used it to stabiUze a WGM resonator at the level of a few 
nano-Kelvin, which will allow us to use this resonator as an ultra-stable laser lock reference. 
The demonstrated technique can be used in a variety of other applications requiring high tem- 
perature stability, as well as ultra-sensitive measurements of temperature variations. As a few 
examples of such applications we would like to suggest the thermal stabilization of quartz os- 
cillator, mid- or far-IR sensitive bolometers, precise calorimetric measurements in chemistry, 
and study of optical and mechanical aging effects in various crystalline resonators. 
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